## A generating matrix of a PRS code

from PyM import *

# def PRS_G(w,k):
#     n = order(w)
#     return matrix([geometric_series(w**j,n) for j in range(k)])

w = primitive_root(Zn(13))

G = PRS_G(w,6)

show( G )